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We report on measurements of the inclusive jet production cross section as a function of the 
jet transverse momentum in pp collisions at v'* = 1-96 TeV, using the ki algorithm and a data 
sample corresponding to 1.0 fb~^ collected with the Collider Detector at Fermilab in Run II. The 
measurements are carried out in five different jet rapidity regions with \y^''*'\ < 2.1 and transverse mo- 
mentum in the range 54 < p'^ < 700 GeV/c. Ncxt-to-lcading order perturbative QCD predictions 
are in good agreement with the measured cross sections. 
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I. INTRODUCTION 

The mcasiircmont of the inclusive jet cross section as a function of the jet transverse momentum, plf^ , in pp colUsions 
at ^/s = 1.96 TeV constitutes a test of perturbative quantum chromodynamics (pQCD) [1]. In Rim II of the Tevatron, 
measurements of the jet cross section for jets with p'^^ up to about 700 GeV/c [2, 3] have extended the p*^^ range by 
more than 150 GeV/c compared to Run I [4-6]. In particular, the CDF collaboration recently published results [2] 
on inclusive jet production using the fcx algorithm [7, 8] for jets with jf^ > 54 GeV/c and rapidity [9] in the region 
0.1 < lyj"*! < 0.7, which are well described by next-to-leading order (NLO) pQCD predictions [10]. As discussed 
in [2], the Ut algorithm has been widely used for precise QCD measurements at both e+e^ and e^p colliders, and 
makes possible a well defined comparison to the theoretical predictions [8]. The pQCD calculations involve matrix 
elements, describing the hard interaction between partons, convoluted with parton density functions (PDFs) [11, 12] 
in the proton and antiproton that require input from experiment. The pQCD predictions are affected by the still- 
limited knowledge of the gluon PDF, which translates into a large uncertainty on the theoretical cross sections at 
high p^^* [2, 3]. Inclusive jet cross section measurements from Run I at the Tevatron [5], performed in different jet 
rapidity regions, have been used to partially constrain the gluon distribution in the proton. This article continues 
the studies on jet production using the /ct algorithm at the Tevatron [2, 6] and presents new measurements of the 
inclusive jet production cross section as a function of ^tJ^* in five different jet rapidity regions up to = 2.1, based 
on 1.0 fb~^ of CDF Run II data. The measurements are corrected to the hadron level [13] and compared to NLO 
pQCD predictions. 

II. EXPERIMENTAL SETUP 

The CDF II detector (see Fig. 1) is described in detail in [14]. The sub-detectors most relevant for this analysis 
are discussed briefly here. The detector has a charged particle tracking system immersed in a 1.4 T magnetic field. 
A silicon microstrip detector [15] provides tracking over the radial range 1.35 to 28 cm and covers the pseudorapidity 
range < 2. A 3.1-m-long open-cell drift chamber [16] covers the radial range from 44 to 132 cm and provides 
tracking coverage for jryj < 1. Segmented sampling calorimeters, arranged in a projective tower geometry, surround 
the tracking system and measure the energy of interacting particles for \'rj\ < 3.6. The central barrel calorimeter [17] 
covers the region \r]\ < 1. It consists of two sections, an electromagnetic calorimeter (CEM) and a hadronic calorimeter 
(CHA), divided into 480 towers of size 0.1 in r] and 15° in (f). The end-wall hadronic calorimeter (WHA) [18] is behind 
the central barrel calorimeter in the region 0.6 < |r/| < 1.0, providing forward coverage out to |7;| < 1.3. In Rim 11, new 
forward scintillator-plate calorimeters [19] replaced the Run I gas calorimeter system. The new plug electromagnetic 
calorimeter (PEM) covers the region 1.1 < |ry| < 3.6, while the new hadronic calorimeter (PHA) provides coverage 
in the 1.3 < \ri\ < 3.6 region. The calorimeter has gaps at jryj « (between the two halves of the central barrel 
calorimeter) and at |?7| « 1.1 (in the region between the WHA and the plug calorimeters). The measured energy 
resolutions for electrons in the electromagnetic calorimeters [17, 19] are 14%/ \/Et 2% (CEM) and 16%/ VE ® 1% 
(PEM), where the energies are expressed in GeV. The singlc-pion energy resolutions in the hadronic calorimeters, as 
determined in test-beam data [17-19], are 50%/^/E^e3% (CHA), 75%/V^®4% (WHA) and 80%//E©5% (PHA). 
Cherenkov counters covering the 3.7 < \r]\ < 4.7 region [20] measure the average number of inelastic pp collisions per 
bunch crossing and thereby determine the beam luminosity. 

III. JET RECONSTRUCTION 

The fc'i' algorithm [8] is used to reconstruct jets from the energy depositions in the calorimeter towers in both 
data and Monte Carlo simulated events (see Section VI). For each calorimeter tower, the four-momenta [21] of its 
electromagnetic and hadronic sections are summed to define a physics tower. First, all physics towers with transverse 
momentum above 0.1 GeV/c are considered as protojets. The quantities 

^T,i = PT,i ; fcT,(i,j) = min(p|^i,PTj) • Ai??j/£)^ (1) 

are computed for each protojet and pair of protojets, respectively, where pT,i denotes the transverse momentum 
of the i**^ protojet, Ai?i,j is the distance {y — (j) space) between each pair of protojets, and D is a parameter that 
approximately controls the size of the jet by limiting, in each iteration, the clustering of protojets according to their 
spacial separation. All fcx.i and A;T,(i,j) values are then collected into a single sorted list. In this list, if the smallest 
quantity is of the type fcx,!, the corresponding protojet is promoted to be a jet and removed from the list. Otherwise, if 
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the smallest quantity is of the type A;T,(i,j), the protojets are combined into a single protojet by summing up their four- 
vector components. The procedure is iterated over protojets until the list is empty. The jet transverse momentum, 
rapidity, and azimuthal angle arc denoted as J^x^^ai' yiai' "^cai' respectively. 

In the Monte Carlo event samples, the same jet algorithm is also applied to the final-state particles, considering all 
particles as protojets, to search for jets at the hadron level. The resulting hadron-level jet variables are denoted as 
<ad, and Cad- 

IV. EVENT SELECTION 

Events are selected online using a three-level trigger system [22] with unique sets of selection criteria called paths. 

For the different trigger paths used in this measurement, this selection is based on the measured energy deposits in 
the calorimeter towers, with different thresholds on the jet Et and different prescale factors [23] (see Table I). In 
the first-level trigger, a single trigger tower [24] with Et above 5 GeV or 10 GeV, depending on the trigger path, is 
required. In the second-level trigger, calorimeter clusters are formed around the selected trigger towers. The events 
are required to have at least one second-level trigger cluster with Et above a given threshold, which varies between 
15 and 90 GeV for the different trigger paths. In the third-level trigger, jets arc reconstructed using the CDF Run I 
cone algorithm [25], and the events are required to have at least one jet with Et above 20 to 100 GeV. 



Trigger Path 


Level 1 tower Ex [GeV] 


Level 2 cluster Et [GeV] 


Level 3 jet Et [GeV] 


eff. prescale 


JET 20 


5 


15 


20 


775 


JET 50 


5 


40 


50 


34 


JET 70 


10 


60 


70 


8 


JET 100 


10 


90 


100 


1 



TABLE I: Summary of trigger paths, trigger thresholds and effective prescale factors employed to collect the data. 

Jets are then reconstructed using the fcx algorithm, as explained in Section III, with D = 0.7. For each trigger path, 
the minimum y^x'cai' each ll/cail region, is chosen in such a way that the trigger selection is fully efficient. The 
efficiency for a given trigger path is obtained using events from a different trigger path with lower transverse energy 
thresholds (see Table I). In the case of the JET 20 trigger path, the trigger efficiency is extracted from additional 
control samples, which include a sample with only first-level trigger requirements as well as data collected using 
unbiased trigger paths with no requirement on the energy deposits in the calorimeter towers. As an example, for jets 
in the region 0.1 < < 0.7, Fig. 2 shows the trigger efficiency as a function ofpl^^^^^ for the different samples. The 
following selection criteria have been imposed: 

1. Events are required to have at least one reconstructed primary vertex with 2;-position within 60 cm of the 
nominal interaction point. This partially removes beam-related backgrounds and ensures a well-understood 

event-by-event jet kinematics. 

2. Events are required to have at least one jet with rapidity in the region ^ -^-1 ^nd corrected Px*cai (^'-'^ 
Section IX) above 54 GeV/c, which constitutes the minimum jet transverse momentum considered in the analysis. 
The measurements are limited to jets with < 2.1 to avoid contributions from the p and p remnants that 
would affect the measured p^T^ai most forward region of the calorimeter. 

3. In order to remove beam-related backgrounds and cosmic rays, the events are required to fulffil ^t/V^Et < 
-^(P'r'cai)' where ^t denotes the missing transverse energy [26] and T,Et = -^t total transverse energy 
of the event, as measured using calorimeter towers with E^ above 0.1 GeV. The threshold function -F(PT°cai) 
defined as F{pi^*^) ~ TOm(2-|-0.0125 xpi^'', 7), where PT^^ai uncorrected transverse momentum of the leading 
jet (highest pi^*) in GeV/c, and F is in GeV^^^. This criterion preserves more than 95% of the QCD events, as 
determined from IVIonte Carlo studies (see Section VI). A visual scan of events with pi^^g^i > 400 GeV/c did not 
show remaining backgrounds. 

[Measurements are carried out in five different jet rapidity regions: < 0.1, 0.1 < |?ycail < ^.7, 0.7 < < l-l, 

1.1 < < 1.6, and 1.6 < ly^^j] < 2.1, where the different boundaries are chosen to reduce systematic effects coming 
from the layout of the calorimeter system. 
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V. EFFECT OF MULTIPLE PP INTERACTIONS 

The measured pJ^*^.^! includes contributions from multiple pp interactions per bunch crossing at high instantaneous 
luminosity. The data used in this measurement wore collected at between 0.2 x lO'^^cm^^s"'^ and 16.3 x 

10'^^cm~^s^^ with an average of 4.1 x lO'^^cm^^s^^. On average, 1.5 inelastic pp interactions per bunch crossing are 
expected. At the highest considered, an average of 5.9 interactions per bunch crossing are produced. This mainly 
affects the measured cross section at low pJ^* where the contributions are sizeable. Multiple interactions are identified 
via the presence of additional primary vertices reconstructed from charged particles. The measured jet transverse 
momenta are corrected for this effect by removing a certain amount of transverse momentum, x {Ny — 1), where 
Ny denotes the number of reconstructed primary vertices in the event and (5™ is determined from the data by 
requiring that, after the correction is applied, the ratio of cross sections at low and high £'"^'' does not show any pi^*' 
dependence. The study is carried out separately in each region, and the results are consistent with a common 
value (5™ = 1.86 ± 0.23 GeV/c across the whole rapidity range. 

VI. MONTE CARLO SIMULATION 

Monte Carlo simulated event samples are used to determine the response of the detector and the correction factors to 
the hadron level. The generated samples are passed through a full CDF II detector simulation (based on GEANtS [27], 
where the GFLASH [28] package is used to simulate the energy deposition in the calorimeters) and then reconstructed 
and analyzed using the same analysis chain as used for the data. 

Samples of simulated inclusive jet events have been generated with PYTHIA 6.203 [29] and HERWIG 6.4 [30] Monte 
Carlo generators, using CTEQ5L [31] PDFs. The PYTHIA samples have been created using a specially tuned set of 
parameters, denoted as PYTHIA-TUNE A [32], that includes enhanced contributions from initial-state ghion radiation 
and secondary parton interactions between remnants. The parameters were determined from dedicated studies of 
the underlying event using the CDF Run I data [33] and has been shown to properly describe the measured jet 
shapes in Run II [34]. In the case of pythia, fragmentation into hadrons is carried out using the string model [35] as 
implemented in jetset [36], while herwig implements the cluster model [37]. 

VII. SIMULATION OF THE CALORIMETER RESPONSE TO JETS 

Dedicated studies have been performed to validate the Monte Carlo simulation of the calorimeter response to jets for 
the different [y^^^*] regions. Previous analyses [2] for jets with 0.1 < < 0.7 indicate that the simulation properly 
reproduces both the average p^^* and the jet momentum resolution, tr te* , as measured in the data. The study is 

performed for the rest of the regions using jets in the range 0.1 < [y^^*] < 0.7 as a reference. An exclusive dijet 
sample is selected, in data and simulated events, with the following criteria: 

1. Events are required to have one and only one reconstructed primary vertex with ^-position within 60 cm of the 

nominal interaction point. 

2. Events are required to have exactly two jets with ^x^cai ^ GeV/c, where one of the jets must be in the region 
0.1 < < 0.7. 

3. ^t/VS-Et < -F(Pt!c"i)' as explained in Section IV. 

The bisector method [38] is applied to data and simulated exclusive dijet events to test the accuracy of the simulated 
cr Jet in the detector. The study indicates that the simulation systematically underestimates the measured a jot by 

6% and 10% for jets in the regions 0.7 < |?/^^*| < 1.1 and 1.6 < < 2.1, respectively, with no significant pJ^^j.^! 

dependence. An additional smearing of the reconstructed p!^\gx is applied to the simulated events to account for 
this effect. In the region 1.1 < Ii^^^yI < 1.6, the measured a^at is overestimated by 5% in the simulation. The effect 

on the final result is included via slightly modified unfolding factors (see Section IX). For jets with |?/^^*| < 0.1, the 
simulation properly describes the measured u jet . Figure 3 shows the ratio between u jet in data and simulated events, 

Pt P'y 

^data/p.mc ^ in different [j/^!^*! regions as a function of the average p!^ of the dijet event. After corrections have 
been applied to the simulated events, data and simulation agree. In the region 1.1 < < 1.6, and only for the 
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purpose of presentation, a 5% smearing of the reconstructed p^x^^ai ^ applied to the data to show the resulting good 
agreement with the uncorrected simulated resolution. The relative difference between data and simulated resolutions 
is conservatively taken to be ±8% (see Fig. 3) over the whole range in pi^^^^i and in the evaluation of systematic 
uncertainties. 

The average jet momentum calorimeter response in the simulation is then tested by comparing the pi^^g;^ balance 
in data and simulated exclusive dijet events. The variable /?, defined as [39] 

, . , test jet ref . jet 
1 + (A) PT,cal ~ PT,cal 
^ = ^ TaT ' ^ = test iet ;^lit ' (2) 

1-(A) PT%i +Pt'cA 

is computed in data and simulated events in bins of (Pxcai°* + ^'t cai''*)/^' where Pxcai*^* denotes the transverse 
momentum of the jet in the region 0.1 < < 0.7, and clr* transverse momentum of the jet in the 

region under study. Figure 4 presents the ratios /3data//3mc as a function of p!^\^i = Pt clr* ™ different bins. 

The study indicates that small corrections arc required around calorimeter gaps, < 0.1 and 1.1 < \l/^ai\ ^ l-^' ^ 
well as in the most forward region, 1.6 < < 2.1. For jets with > 1.1, the correction shows a moderate pi^^^^i 

dependence, and several parameterizations are considered to extrapolate to very high pi^^^^i- The difference observed 
in the final results, using different parameterizations, is included as part of the total systematic uncertainty. 



VIII. RECONSTRUCTION OF THE JET VARIABLES 



The jet reconstruction in the detector is studied using Monte Carlo event samples, with modified jet energy response 
in the calorimeter, as described in the previous section, and pairs of jets at the calorimeter and hadron levels matched 

in (y - (f)) space by requiring (j/^^* — yjf^^)^ + {(j^^cai ~ '^had)^ < ^- These studies indicate that the angular variables 
of a jet are reconstructed with no significant systematic shift and with a resolution better than 0.05 units in y and 
(/> at low Px*cai' improving as Px*cai increases. The measured Px*cai systematically underestimates that of the hadron 
level jet. This is attributed mainly to the non-compensating nature of the calorimeter [17]. For jets with p^x*^,^! around 
50 GeV/c, the jet transverse momentum is reconstructed with an average shift that varies between —9% and —30% 
and a resolution between 10% and 16%, depending on the region. The jet reconstruction improves as Pxcai 

increases. For jets with p^^^^^^ around 500 GeV/c, the average shift is —7% and the resolution is about 7%. 



IX. UNFOLDING 

The measured pi^^^^^ distributions in the different regions are unfolded back to the hadron level using simulated 
event samples (see Section VI), after including the modified jet energy response described in Section VII. pythia- 
TUNE A provides a reasonable description of the different jet and underlying event quantities, and is used to determine 
the correction factors in the unfolding procedure. In order to avoid any potential bias on the correction factors due 
to the particular PDF set used during the generation of the simulated samples, which translates into slightly different 
simulated p^^^^ai distributions, the underlying pt spectrum [40] in PYTHIA-TUNE A is re- weighted until the Monte Carlo 
samples acciirately follow each of the measured Px*cai distributions. The unfolding is carried out in two steps. 

First, an average correction is computed separately in each jet rapidity region using corresponding matched pairs 
of jets at the calorimeter and hadron levels. The correlation (p^^^^d ~^'T*cai) versus (Pr^cai) (^^^ ^^S- 5), computed in 
bins of (Pt had +PT*cai)/^' uscd to extract correction factors which are then applied to the measured jets to obtain 
the corrected transverse momenta, Pxrar- ^ eSiCh. jet rapidity region, a cross section is defined as 

d'a _ 1 TVjf, 

where TVjg* denotes the number of jets in a given p^^^^-oi. bin, Ap!^^^^^ is the size of the bin, Ay^^* denotes the size of 
the region in t/^^*, and £ is the integrated luminosity. N^^^ includes event-by-event weights that account for trigger 
prescale factors, and ApJ^^^^^. is chosen according to the jet momentum resolution. 
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Second, each measurement is corrected for acceptance and smearing effects using a bin-by-bin unfolding procedure, 
which also accounts for the efficiency of the selection criteria. The unfolding factors, defined as 



c^(pv"co.c^) = ./J r;r . (4) 



are extracted from Monte Carlo event samples and applied to the measured p!^ distributions to obtain the final 
results. As shown in Fig. 6, the factor t^(PT*cor' ^03*1) increases with p!^^^^^ and presents a moderate dependence. 
At low p^x*cor' unfolding factor varies between 1.02 and 1.06 for different rapidity regions. For jets with Px*cor 
of about 300 GeV/c, the factor varies between 1.1 and 1.2, and increases up to 1.3 - 1.4 at very high jJ^^o^- 
region 1.1 < < 1.6, the unfolding factor includes an additional correction, /u(PT*cor)' *o account for the fact 
that the simulation overestimates the jet momentum resolution in that region (see Section VII). The factor fviPT^cor) 
is computed from Monte Carlo samples as the ratio between the p^x*had distribution smeared using the simulated a^jat 
and the one smeared using cr^ct in data as extracted from the bisector method (see Section VII). The factor /u(f'T*cor) 
is about 1.03 and shows no significant j3^x*cor dependence. 



X. SYSTEMATIC UNCERTAINTIES 



A detailed study of the systematic imccrtainties on the measurements has been carried out [41]. Tables II-III show 
the different contributions to the total systematic uncertainty in each pl^^ bin and ly-'"*! region: 

1. The measured jet energies are varied by ±2% at low pi^* to ±2.7% at high pl^^ to account for the uncertainty on 
the absolute energy scale in the calorimeter [42] . This introduces an uncertainty on the measured cross sections 
which varies between ±9% at low pj^* and ^39% at high pl^*, and dominates the total systematic uncertainty on 
the different measurements. 

2. Several sources of systematic uncertainty on the ratio /3data//3inc ai"© considered for the different jj/-'''*] regions: 

(a) The uncertainty on the definition of the exclusive dijet sample in data and Monte Carlo events introduces 
a ±0.5% uncertainty on the absolute energy scale for jets outside the region 0.1 < < 0.7, which 
translates into an uncertainty on the cross sections between ±2% at low p!^^ and ±10% at very high p!^^. 

(b) The use of different /3data//3mc parameterizations for jets with \y'^^\ > 1.1 introduces uncertainties between 

12% and 23% at very high 

(c) In the region 1.1 < < 1.6, an additional uncertainty on the measured cross sections, independent 

of p!^^, accounts for variations in the /3data//3mc ratio due to the overestimation of the jet momentum 
resolution in the simulated samples. 

3. A ±8% uncertainty on the jet momentum resolution introduces an uncertainty between ±2% at low p!^*' and 

athighpjp^*. 

4. The unfolding procedure is repeated using herwig instead of pythia-tune a to accoimt for the imccrtainty on 
the modeling of the parton cascades and the jet fragmentation into hadrons. This translates into an uncertainty 
on the measured cross sections between ±3% and ±8% at low pJ^* that becomes negligible at very high pi^*. 

5. The unfolding procedure is also carried out using unweighted pythia-tune a, to estimate the residual depen- 
dence on the p^r^* spectra. This introduces an uncertainty of about ±3% to ±7% at very high pJ^*, which becomes 

negligible at low p!^* . 

6. The quoted ±0.23 GeV/c uncertainty on (5™ is taken into account. The maximal effect on the measured cross 
sections is about ±2%. 

7. Different sources of systematic uncertainty related to the selection criteria are considered. The threshold on the 
2-position of the primary vertex is varied by ±5 cm in data and simulated events. The lower edge of each pi^^^^i 
bin is varied by ±3% in data and simulated events. The scale is varied by ±10% in the data. The total 
effect on the measured cross sections is smaller than 1% and considered negligible. 
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Positive and negative deviations with respect to the nominal values in each p!^* bin are added separately in quadrature. 
Figure 7 shows the total systematic uncertainty as a function of pJ^* in the different \y-''^^\ regions, where an additional 
5.8% uncertainty on the total luminosity is not included. 

XI. QCD PREDICTIONS 

The measurements are compared to parton-level NLO pQCD predictions, as computed using jetrad [10] with 
CTEQ6.1M PDFs [11] and the renormaHzation and factorization scales {^r and iip) both set to = max{p!^^)/2. 
Different sources of uncertainty on the theoretical predictions are considered. The main contribution comes from 
the uncertainty on the PDFs and is computed using the Hessian method [43]. At low pJ^*, the uncertainty is about 
±5% and approximately independent of 7/^°*. The uncertainty increases as p!^^ and increase. At very high pJ^*, 
the uncertainty varies between a-nd ^40%" for jets with \y-'°*'\ < 0.1 and 1.6 < ly-'®*] < 2.1, respectively, and 

is dominated by the limited knowledge of the gluon PDF. An increase of iir and /j-f from iJ.o to 2iJ.o changes the 
theoretical predictions by only a few percent. Values significantly smaller than /xq lead to unstable NLO results and 
are not considered. 

The theoretical predictions include a correction factor, Chad(?5t*' 2/^*^*)' approximately accounts for non- 

perturbative contributions from the underlying event and fragmentation of partons into hadrons (sec Fig. 8 and 
Tables IV- V). In each jet rapidity region, Chad is estimated, using PYTHIA-tune a, as the ratio between the nominal 
Pt. ,had distribution and the one obtained after removing the interactions between p and p remnants and the fragmen- 
tation into hadrons in the Monte Carlo samples. The correction decreases as pJ^* increases and shows a moderate 
Ij/j*^*! dependence. At low pi^*, Chad varies between 1.18 and 1.13 as |t/j''*| increases, and it becomes of the order of 
1.02 at very high pl^^. The uncertainty on Chad varies between ±9% and ±12% at low pJ^* and decreases to about 
±1% at very high p^^*, as determined from the difference between the parton-to-hadron correction factors obtained 
using HERWIG instead of pythia-tune a. 

XII. RESULTS 

The measured inclusive jet cross sections, (Pa/dp!^^dy^''*', refer to hadron-level jets, reconstructed using the fcx 
algorithm with D — 0.7, in the region pi^* > 54 GeV/c and \y^'='^\ < 2.1. Fi gure 9 shows the measured cross sections 
as a function of pi^* in five different regions compared to NLO pQCD predictions. The data are reported in 

Tables IV-V. The measured cross sections decrease by more than seven to eight orders of magnitude as pi^^ increases. 
Figure 10 shows the ratio data/theory as a function of pJ^* in the five different \y^'^^\ regions. Good agreement is 
observed in the whole range in p^^' and yj®* between the measured cross sections and the theoretical predictions. 
In particular, no significant deviation from the pQCD prediction is observed for central jets at high p^^'. The 
corresponding tests, relative to the nominal pQCD prediction and performed separately in each region, give 
probabilities that vary between 9% and 90%. A global test, applied to all the data points in all jt/j*^*! regions 
simultaneously, gives a probability of 7%. In both cases, a detailed treatment of correlations between systematic 
uncertainties was considered, as discussed in Appendix A. In addition. Fig. 10 shows the ratio of pQCD predictions 
using MRST2004 [12] and CTEQ6.1M PDF sets, well inside the theoretical and experimental uncertainties. The 
uncertainty on the measured cross sections at high p^^*, compared to that on the theoretical predictions, indicates 
that the data presented in this article will contribute to a better understanding of the gluon PDF. 

Finally, in the region 0.1 < < 0.7, the analysis is repeated using different values for D in the fcx algorithm: 

D = 0.5 and D = 1.0. In both cases, good agreement is observed between the measured cross sections and the 
NLO pQCD predictions in the whole range in pJ^* (see Fig. 11 and Tables VI- VII). The corresponding tests give 
probabilities of 84% and 22% for D = 0.5 and D = 1.0, respectively. As D decreases, the measurement is less 
sensitive to contributions from multiple pp interactions per bunch crossing, and the presence and proper modeling of 
the underlying event. For D = 0.5 {D = 1.0), the value for 5™ becomes 1.18 ± 0.12 (3.31 ± 0.47) GeV/c, and the 

parton-to-hadron correction factor applied to the pQCD predictions is Chad = 1.1 (Chad = 1-4) at low pJ^*. 
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XIII. SUMMARY AND CONCLUSIONS 



We have presented results on inclusive jet production in pp collisions at y/s = 1.96 TeV for jets with transverse 

momentum pl^- > 54 GoV/c and rapidity in the region < 2.1, using the fcx algorithm and based on 1.0 fb^^ of 

CDF Run II data. The measured cross sections are in agreement with NLO pQCD predictions after the necessary 
non-perturbative parton-to-hadron corrections are taken into account. The results reported in this article should 
contribute to a better understanding of the gluon PDF inside the proton when used in QCD global fits. 
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FIG. 2: Measured trigger efficiencies as a function of p^x*cai ^'^^ different trigger paths and in the region 0.1 < < 0.7. 
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FIG. 3: Ratio cr'^jvt''/cr™5t as a function of the average Px'cai of ^'^'^ dijet event, in different \y'^^i\ regions, before (black squares) 

and after (open circles) corrections have been applied (see Section VII). The solid lines are fits to the corrected ratios. The 
dashed lines indicate a ±8% relative variation considered in the study of systematic uncertainties. 




FIG. 4: Ratio /3data//3mc as a function of pi^'^^^j in different regions. The solid lines show the nominal parameterizations 

based on fits to the ratios. In the region > 1.1, the dashed lines indicate different parameterizations used to describe the 
ratios at high pi^^^^, and are considered in the study of systematic uncertainties. 



14 












r iyiALi<o.i 




0.1<|yi^;j<0.7 












: ' ' ' ' 


: ■ ■ ■ " ■ 


- 1 . . 


' 1 . . 


: 0.7<|yi^;i<1.1 




: 1.1<|yi^;j<l.6 _ 














■ , , , 1 , , , 1 , , , 1 , , 


\ 1-6<|yi^J<2.1 




200 


400 600 






PrSoplGeV/c] 


















- , , , 1 , , , 1 , , , 1 , , 







200 400 600 



FIG. 6: Unfolding factors, f^(p^T*cori fiS)! extracted from PYTHIA TUNE A simulated event samples, as a function of pl^^^^ in 
the different 12/^2 1 regions. 
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FIG. 7: Total systematic uncertainty (in percent) on the measured inclusive differential jet cross sections as a function pi^* for 
the different It/^^l regions (see Tables II-III). An additional 5.8% uncertainty on the integrated luminosity is not included. 
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FIG. 8: Magnitude of the parton-to-hadron correction, CHAD(pi^', J/^"^*), used to correct the NLO pQCD predictions (see 
Tables IV-V). The shaded bands indicate the quoted Monte Carlo modeling uncertainty. 
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TABLE II: Systematic uncertainties (in percent) on the measured inclusive jet differential cross section as a function of pJ^* for 
jets in the regions ly'^'^l < 0.1 and 0.1 < ly'"*] < 0.7 (see Fig. 7). The different columns follow the discussion in Section X. An 
additional 5.8% uncertainty on the integrated luminosity is not included. 
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TABLE III: Systematic uncertainties (in percent) on the measured inclusive jet differential cross section as a function of 
for jets in the regions 0.7 < < 1.1, 1.1 < < 1.6, and 1.6 < < 2.1 (see Fig. 7). The different columns follow the 
discussion in Section X. An additional 5.8% uncertainty on the integrated luminosity is not included. 
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X 10" 
X 10" 
X 10" 



1.177 ±0.124 
1.144 ± 0.097 
1.119 ±0.077 
1.098 ± 0.061 
1.083 ±0.049 
1.070 ± 0.039 
1.060 ± 0.032 
1.052 ±0.026 
1.046 ±0.021 
1.041 ±0.017 
1.037 ±0.013 
1.034 ±0.010 
1.032 ± 0.007 
1.030 ±0.005 
1.028 ±0.002 
1.027 ± 0.001 
1.026 ±0.006 



(0.1 < < 0.7) 



ict 

[GeV/c] 



a ± (stat.) 



± (sys. 
GeY/c] 



Chad 
parton — > hadron 



54 - 62 

62 - 72 
72 - 83 
83 - 96 
96 - 110 
110 - 127 
127 - 146 
146 - 169 
169 - 195 
195 - 224 
224 - 259 
259 - 298 
298 - 344 
344 - 396 
396 - 457 
457 - 527 
527 - 700 



(14.0±0.20lJ;° 

6.14 ±0.12 
2.69 ± 0.02 



+0.66 
-0.65 
0.29 
0.27 

l.l4±0.0ll°;J^ 

4.90 ± 0.04to |s 

2.08 ± 0.02to;2i 
8.51±0.04l°;8B 

3.33 ± 0.0210-4? 
1.23 ±0.01+0-16 

4.53 ± 0.02l°;6f 

+0.26 
0.24 

4.87 ±0.0610-91 

1.43 ± 0.02t° 31 

3.69±0.10+°-|4 

7.18±0.34+j-g0 
1.16±0.13l°;44 
8.97 ± 2.40l^;^^ 



1.57 ±0.01^ 



X 10" 

X loO 
X loO 
X 10° 
X 10" 
X 10" 
X 10" 
X 10" 
X 10- 
X 10" 
X 10- 
X 10" 
X 10" 
X 10" 
X 10" 
X 10" 
X 10" 



1.188 ±0.140 

1.156 ±0.113 
1.129 ±0.091 
1.108 ± 0.073 
1.090 ±0.059 
1.076 ±0.047 
1.065 ± 0.038 
1.055 ±0.029 
1.047 ±0.023 
1.041 ±0.017 
1.036 ±0.012 
1.031 ±0.007 
1.028 ± 0.003 
1.025 ±0.001 
1.023 ± 0.004 
1.021 ± 0.008 
1.018 ±0.014 



TABLE IV: Measured inclusive jet differential cross section as a function of pl^* for jets in the regions \y^''*^\ < 0.1 and 
0.1 < < 0.7 (see Fig. 9). An additional 5.8% uncertainty on the integrated luminosity is not included. The parton-to- 
hadron correction factors, Chad(p^x*, 2/-'''*), are applied to the pQCD predictions (see Fig. 8). 
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(0.7 < |2/j<=*| < 1.1) 



[GeV/c] 



a ± (stat.) ± (sys.) 
[nb/(GeV/c)] 



Chad 
parton — » hadron 



54 - 62 
62 - 72 
72 - 83 
83 - 96 
96 - 110 
110 - 127 
127 - 146 
146 - 169 
169 - 195 
195 - 224 
224 - 259 
259 - 298 
298 - 344 
344 - 396 



396 
457 



457 
527 



(I2.3±0.2l^;;;) 
2.40±0.02l°;2^) 

i-ooio-oiif,:;;;;) 

0.46) 

1.73±0.03«Ji) 
6.83 ±0.051°;!^) 

2.52±0.03l°:33) 

3.04 ± omtolj) 

9.52±0.1lt};f|) 
2.53 ±0.05^^;!^) 
6.18±0.17lJ;«^) 
l.ll±0.07l°-j^) 

i..53±o.2ot;;;5-;;) 



2.17 = ().■ 



X 10" 
X 10" 
X 10° 
X 10" 
X 10^1 
X 10-1 
X 10-2 
X 10-2 
X 10-3 
X 10-3 
X 10-* 
X 10-* 
X 10-5 
X 10-5 
X 10-'' 
X 10 " 



1.169 ±0.125 
1.143 ±0.103 
1.120 ± 0.085 
1.102 ±0.070 
1.087 ±0.057 
1.075 ±0.047 
1.064 ± 0.038 
1.056 ±0.031 
1.048 ±0.024 
1.042 ± 0.019 
1.037 ±0.014 
1.033 ±0.009 
1.030 ± 0.005 
1.027 ±0.001 
1.025 ± 0.003 
1.02:-! = 0.007 



(1.1 < < 1.6) 



[GeV/c] 


54 


- 62 


62 


- 72 


72 


- 83 


83 


- 96 


96 


110 


110 


- 127 


127 


- 146 


146 


- 169 


169 


- 195 


195 


- 224 


224 


- 259 


259 


- 298 


298 


- 344 


311 


- :-!!)() 



a ± (stat.) ± (sys.) 

[nb/(GeV/c)] 



Chad 
parton — » hadron 



(11.0 ± 
(4.40 ± 0. 

(1.82 ± 0. 
(7.22 ±0. 
(2.98 ± 0. 
(1.14 ±0. 
(4.10 ±0. 
(1.39 ±0 
(4.19 ±0 
(1.15 ±0 
(2.73 ± 
(5.18 ±0 
(7.99 ± 0. 
(1.05 = 0. 



.02j 
.04j 
.02^ 
.09j 

.23:* 



3111) 


X 


10° 




1 


160 ±0.125 


,+0.64 X 
'-0.53'' 


X 


10° 




1 


133 ±0.101 


,+0.22x 
'-0.22/ 


X 


10" 




1 


111 ±0.081 


,+0.90x 
-0.90.* 


X 


10- 


1 


1 


094 ± 0.065 


,+0.38 -> 
'-0.38/ 


X 


10- 


1 


1 


080 ± 0.052 


^^^) 


X 


10- 


1 


1 


068 ± 0.042 


,+0.60\ 
'-0.60^ 


X 


10- 


2 


1 


059 ± 0.034 


)+0.22x 
'-0.23/ 


X 


10- 


2 


1 


051 ± 0.027 


1+0.78 N 
'-0.76^ 


X 


10- 


3 


1 


045 ± 0.021 


)+0.25x 
'-Q.2A> 


X 


10- 


3 


1 


040 ± 0.016 


.+0.73-, 
'-0.64-' 


X 


10- 


4 


1 


036 ± 0.012 


>+1.68>, 

'-1.39'' 


X 


10- 


5 


1 


033 ± 0.009 


15;S) 


X 


10- 


(i 


1 


030 + 0.006 




X 


10 


1) 


1 


028 = 0.003 



(1.6 < < 2.1) 



[GeV/c] 



(T ± (stat.) ± (sys.) 

[nb/(GeV/c)] 



Chad 

parton hadron 



54 - 62 
62 - 72 
72 - 83 
83 - 96 
96 - 110 
110 - 127 
127 - 146 
146 - 169 
169 - 195 
195 - 224 
224 - 259 
259 - 298 



6.67±0.15l2;**) 
2.68 ± 0.02t°;^2) 
1.04±0.0ll°;12) 
:3.77±0.04«;«) 

1.32±0.02l° 

4.18±0.04t°-g) 

1.21±0.02l0;24) 
2.92 ±0.041°;^?) 

5.74±0.09tl;||) 

8.49 ±0.3112:42) 
8.65±0.63l3;Jg) 
:5.67±1.65l^;2g) 



X 10" 
X 10° 

X 10° 
X 10-1 
X 10-1 

X 10-2 

X 10-2 

X 10-3 
X 10-* 
X 10-5 
X 10-s 
X 10-'^ 



1.132 ±0.104 
1.116 ±0.087 
1.100 ±0.072 
1.086 ±0.058 
1.072 ± 0.045 
1.059 ±0.033 
1.047 ±0.022 
1.035 ±0.012 
1.024 ±0.003 
1.013 ±0.005 
1.003 ± 0.012 
0.993 ± 0.018 



TABLE V: Measured inclusive jet difFcrcntial cross section as a function of pi^' for jets in the regions 0.7 < < 1.1, 

1.1 < Ij/^'^*! < 1.6, and 1.6 < jy^'^'l < 2.1 (see Fig. 9). An additional 5.8% uncertainty on the integrated luminosity is not 
included. The parton-to-hadron correction factors, Chad(p^^*, j/i'^*), are applied to the pQCD predictions (see Fig. 8). 
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Systematic uncertainties [%] (0.1 < < 0.7) {D = 0.5) 



pj^* [GeV/c] 


jet energy scale 


resolution 


unfolding 


pi^* -spectra 


^mi 

Pt 


54 


- 62 


+9.9 
-9.2 


+2.4 
-2.3 


±5.4 


±0 


6 


+ 0.8 
-0.8 


62 


- 72 


+9.8 
-9.0 


+2.4 
-2.2 


±4.8 


±0 


6 


+0.7 
-0.7 


72 


- 83 


+9.8 
-8.9 


+2.3 
-2.2 


±4.3 


±0 


6 


+0.6 
-0.7 


83 


- 96 


+9.7 
-8.9 


+2.2 
-2.1 


±3.8 


±0 


6 


+0.6 
-0.6 


96 


110 


+9.8 
-9.0 


+2.2 
-2.1 


±3.4 


±0 


6 


+0.5 
-0.5 


110 


- 127 


+ 10.0 

-9.4 


+2.1 

-2.0 


±3.1 


±0 


6 


+0.5 
-0.5 


127 


- 146 


+ 10.4 
-9.9 


+2.1 
-2.0 


±2.8 


±0 


6 


+0.4 
-0.4 


146 


- 169 


+ 11.2 
-10.8 


+2.1 
-2.0 


±2.5 


±0 


5 


+0.4 
-0.4 


169 


- 195 


+12.6 
-11.9 


+2.1 
-2.1 


±2.3 


±0 


4 


+0.4 
-0.4 


195 


- 224 


+14.3 
-13.3 


+2.2 
-2.2 


±2.1 


±0 


3 


+0.4 
-0.3 


224 


- 259 


+16.6 
-16.0 


+2.4 
-2.4 


±1.9 


±0 


2 


+ 0.3 
-0.3 


259 


- 298 


+19.3 
-17.0 


+2.7 
-2.7 


±1.8 


±0 


1 


+0.3 
-0.3 


298 


- 344 


+22.3 
-19.4 


+3.1 
-3.2 


±1.6 


±0 


1 


+0.3 
-0.3 


344 


- 396 


+25.7 
-22.1 


+3.7 
-3.8 


±1.5 


±0 


2 


+0.3 
-0.3 


396 


- 457 


+30.7 
-26.5 


+4.5 
-4.6 


±1.4 


±0 


5 


+0.3 
-0.3 


457 


- 527 


+39.6 
-29.7 


+6.5 
-5.6 


±1.3 


±1 


3 


+0.3 
-0.2 


527 


- 700 


+52.6 
-37.7 


+7.4 

-7.3 


±1.2 


±4 


2 


+0.3 
-0.2 



Systematic uncertainties [%] (0.1 < \yi''^\ < 0.7) (D = 1.0) 



-j^* [GeV/c] 


jet energy scale 


resolution 


unfolding 


pij;* -spectra 


<5™' 

Pt 


54 - 


62 


+10.7 
-9.4 


+2.7 
-2.7 


±5 


6 


±0 


4 


+ 3.5 
-2.9 


62 - 


72 


+10.4 
-9.3 


+2.6 
-2.6 


±4 


9 


±0 


4 


+3.0 
-2.6 


72 - 


83 


+10.3 
-9.2 


+2.4 
-2.4 


±4 


2 


±0 


4 


+2.6 
-2.4 


83 - 


96 


+10.2 
-9.2 


+2.3 
-2.3 


±3 


7 


±0 


4 


+2.3 
-2.2 


96 - 


110 


+10.2 
-9.3 


+2.2 
-2.2 


±3 


2 


±0 


4 


+2.1 
-2.0 


110 - 


127 


+10.4 
-9.6 


+2.1 
-2.1 


±2 


8 


±0 


4 


+ 1.9 
-1.8 


127 - 


146 


+10.8 
-10.1 


+2.0 
-2.0 


±2 


5 


±0 


4 


+ 1.7 
-1.7 


146 - 


169 


+ 11.5 
-10.8 


+1.9 
-1.9 


±2 


1 


±0 


4 


+ 1.6 
-1.6 


169 - 


195 


+ 12.6 
-11.8 


+ 1.9 

-2.0 


±1 


9 


±0 


4 


+ 1.5 
-1.4 


195 - 


224 


+ 13.9 
-13.1 


+ 1.9 

-2.0 


±1 


6 


±0 


3 


+ 1.4 
-1.3 


224 - 


259 


+ 15.8 
-14.7 


+2.1 
-2.2 


±1 


4 


±0 


3 


+ 1.3 
-1.3 


259 - 


298 


+ 18.0 
-16.6 


+2.4 
-2.5 


±1 


3 


±0 


2 


+ 1.3 
-1.2 


298 - 


344 


+20.8 
-18.8 


+2.8 
-2.9 


±1 


1 


±0 


2 


+ 1.2 
-1.1 


344 - 


396 


+24.5 
-21.4 


+3.4 

-3.6 


±1 





±0 


2 


+ 1.2 
-1.1 


396 - 


457 


+30.1 
-24.7 


+4.3 
-4.4 


±0 


8 


±0 


5 


+ 1.1 
-1.0 


457- 


527 


+38.8 
-29.5 


+5.4 
-5.4 


±0 


7 


±1 


1 


+1.1 
-1.0 


527 - 


700 


+49.8 
-37.6 


+7.3 
-7.2 


±0 


6 


±3 


4 


+1.0 
-0.9 



TABLE VI: Systematic uncertainties (in percent) on the measured inclusive jet differential cross section as a function of pi^*, 
for jets in the region 0.1 < \y''^^\ < 0.7 and using D = 0.5 and D = 1.0 (see Fig. 11). The different columns follow the discussion 
in Section X. An additional 5.8% uncertainty on the integrated luminosity is not included. 
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dp, 



- (0.1 < W^'\ < 0.7) {D = 0.5) 





a ± (stat.) ± (sys.) 


[GeV/c] 


[nb/( 


GeV/c)] 


54 - 


62 




X 


1 nO 
10 


62 - 


72 


(4.81 ± 0.03+^:^^ 


X 


10° 


72 - 


83 


/"o no _i_ n ni +0.23 
(^z.uy ± U.U1_Q 22 


X 


10° 


83 - 


96 


m _i_ n ni +0.1(3 


X 


10" 


yo - 


1 1 n 
ilU 


(^o.yo it U.U4 Q 


X 


in-i 


110 - 


127 


(1.71 ± 0.02Iq ^ 


X 


10-1 


127 - 


146 


m 71 4- n ni"t"^-^® 
l^u. / i It U.Ui Q 




IQ-i 


146 - 


169 


/o _l_ n no~t~0.32 
(2.76 ± O.UZ q .^-^ 


X 


10-2 


169 - 


195 


(1.04±0.0lll^;i4 


X 


10-2 


195 - 


224 


(3.87 ± 0.02j;°;^^ 


) X 


10-3 


224 - 


259 


(1.34±O.Olto2? 


X 


10-3 


259 - 


298 


(4.26 ± 0.04to;74 


X 


10-4 


298 - 


344 


(1.22 ± 0.02l° 24 


X 


10-4 


344 - 


396 


(3.16±0.09l°-7i 


X 


10-5 


396 - 


457 


(6.30±0.32tJ;^^ 


X 


10-'^ 


457- 


527 


(1.01±0.12l°;4° 


X 


10-s 


527 - 


700 


(0.83 ± 0.23^0 32 


X 


10-7 



Chad 
parton — » hadron 



1.089 ±0.104 
1.076 ± 0.086 
1.064 ± 0.070 
1.055 ± 0.057 
1.047 ±0.047 
1.041 ±0.037 
1.035 ± 0.029 
1.030 ±0.023 
1.026 ±0.017 
1.022 ±0.012 
1.019 ±0.008 
1.017 ±0.005 
1.015 ±0.002 
1.013 ±0.001 
1.011 ±0.002 
1.010 ±0.003 
1.008 ±0.005 



d^a 



(0.1 < < 0.7) {D = 1.0) 



ict 

[GeV/c] 



a ± (stat.) ± {sys.) 
[nb/(GeV/c)] 



Chad 
parton — > hadron 



54 - 62 

62 - 72 
72 - 83 
83 - 96 
96 - 110 
110 - 127 
127 - 146 
146 - 169 
169 - 195 
195 - 224 
224 - 259 
259 - 298 
298 - 344 
344 - 396 
396 - 457 
457 - 527 
527 - 700 



(20.0 ±0.212.3) 

(8.65±0.04+^;J) 
3.59±0.02lg;3g) 

lA9±O.Olt°o.\l) 
6.27±0.05l°;7°) 

2.63±0.03lo;27) 
1.05±0.0ll° ;J2) 

4.04±0.03l°;4^) 
1.48 ±0.01+°;^^) 
5.41±0.02l0;77) 
1.86±0.0lt°-J°) 
5.77±0.04lJ;°^) 
1.70±0.02l°;|^) 
4.26 ±0.10+0 33) 
8.17 ±0.3612-49) 

1.39±0.14l°;||) 
1.19 ±0.271°;^°) 



X 10° 

X 10° 
X 10° 
X 10° 
X 10-1 
X 10-1 
X 10-1 

X 10^2 
X 10^2 
X 10-3 
X 10-3 
X 10-4 
X 10-4 

X 10-5 
X lO-'' 
X 10-8 
X 10-7 



1.372 ± 0.227 

1.296 ±0.171 
1.236 ±0.129 
1.190 ± 0.098 
1.155 ±0.075 
1.127 ±0.057 
1.105 ±0.044 
1.088 ±0.034 
1.075 ±0.026 
1.065 ±0.019 
1.057 ±0.013 
1.050 ± 0.008 
1.045 ± 0.003 
1.041 ±0.003 
1.038 ± 0.009 
1.036 ±0.015 
1.033 ± 0.027 



TABLE VII: Measured inclusive jet differential cross section as a function of pi^* for jets in the region 0.1 < ly^"^] < 0.7 
using D = 0.5 and D = 1.0 (see Fig. 11). An additional 5.8% uncertainty on the integrated luminosity is not included. The 
parton-to-hadron correction factors, Chad(23J^*), are applied to the pQCD predictions. 
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APPENDIX A: CORRELATIONS OF SYSTEMATIC UNCERTAINTIES 

The correlations among systematic uncertainties in different pl^^ bins and \y^^^\ regions are studied in detail. The 
uncertainty on the absolute jet energy scale is decomposed into different sources considered independent but fully 
correlated across p!^^ bins and \y^'^^\ regions. A ±1.8% uncertainty on the absolute energy scale, independent of p^^*, 
results from the sum in quadrature of four different contributions [42]: a ±0.5% uncertainty from the calorimeter 
stability versus time, a ±1.0% uncertainty due to the modeling of the jet fragmentation, a ±0.5% uncertainty from the 
simulation of the electromagnetic calorimeter response, and a ±1.3% uncertainty from the simulation of the calorimeter 
response at the boundary between calorimeter towers. Other contributions to the absolute energy scale uncertainty 
come from the description of the calorimeter response to hadrons for different ranges in hadron momentum [42]. 
Table VIII shows the resulting relative contributions to the quoted systematic uncertainty on the measured cross 
sections related to the absolute jet energy scale uncertainty. 



pip'* 


pj^* independent 




response to hadrons 




[GeV/c] 


uncertainty 


p < 12 GeV/c 


12 < p < 20 GeV/c 


p > 20 GeV/c 


54 - 62 


90.3 


37.8 


15.2 


13.5 


62 - 72 


90.2 


35.2 


16.1 


19.1 


72 - 83 


89.9 


31.9 


17.0 


24.6 


83 - 96 


89.2 


28.8 


17.3 


30.1 


96 - 110 


88.0 


26.0 


16.9 


35.8 


110 - 127 


86.4 


22.7 


16.4 


41.9 


127 - 146 


84.3 


20.0 


15.1 


47.7 


146 - 169 


82.1 


17.2 


14.1 


52.6 


169 - 195 


79.8 


14.6 


12.7 


57.0 


195 - 224 


77.6 


12.5 


11.5 


60.7 


224 - 259 


75.7 


10.7 


10.3 


63.6 


259 - 298 


73.8 


9.1 


9.2 


66.2 


298 - 344 


72.1 


7.8 


8.2 


68.3 


344 - 396 


70.5 


6.8 


7.3 


70.2 


396 - 457 


69.2 


5.8 


6.4 


71.7 


457 - 527 


68.0 


5.0 


5.7 


72.9 


527 - 700 


66.8 


4.2 


5.0 


74.2 



TABLE VIII: Relative contributions (in percent) to the quoted systematic uncertainty on the measured cross sections related 
to the absolute jet energy scale uncertainty. The second column corresponds to a ±1.8% uncertainty on the absolute energy 
scale, £is discussed in the text. Sources are considered independent and fully correlated in pi^' and \y^'^^\- 



The rest of the systematic uncertainties on the measured cross sections, including that on the total integrated lumi- 
nosity, are also assumed to be independent and fully correlated across pl^^ bins and [yj®*! regions, except those related 
to the /3data//3mc ratio, for which uncertainties in different [y-'''*! regions are uncorrelated. 
A global test is performed according to the formula 



is)? 



^ [Serf]-' + [Jaf (s)]2 



Si 



(Al) 



where is the measured cross section for data point j, crj'*(s) is the corresponding prediction, and s denotes the 
vector of standard deviations, Sj, for the different independent sources of systematic uncertainty. The values for 

cr*'*(s) arc obtained from the nominal NLO pQCD prediction, where s includes the uncertainty on Chad but does 

not consider PDF uncertainties. The uncertainty on Chad is assumed to be fully correlated across pi^* bins and [j/j"*! 
regions. The sums in Eq. (Al) run over 76 data points and 17 independent sources of systematic uncertainty, and 

the is minimized with respect to s. Correlations among systematic uncertainties arc taken into account in (Tj^{s). 
As an example, for a given source i, variations of Si will coherently affect all the (7j^{s) values if the corresponding 
systematic uncertainties are considered fully correlated across p!^^ bins and \y^^^\ regions. 



